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\ Abstract 

We study deep inelastic scattering structure functions from hadrons using different holo- 
graphic dual models which describe the strongly coupled regime of gauge theories in the 
large N limit. Particularly, we consider scalar and vector mesons obtained from holographic 
descriptions with fundamental degrees of freedom, corresponding to Af = 2 supersymmetric 
and non-supersymmetric Yang-Mills theories. We explicitly obtain analytic expressions for 
« ■ the full set of eight structure functions, i.e., Fx, F 2 , g±, g 2 , &i, b 2 , b 3 , 64, arising from the 
£sj ■ standard decomposition of the hadronic tensor of spin-one hadrons. We obtain the relations 
2F% = F 2 and 2b\ = b 2 . In addition, we find 61 ~ 0(F\) as suggested by Hoodbhoy, Jaffe 
r—i . and Manohar for vector mesons. Also, we find new relations among some of these structure 
j> ! functions. 



1 Introduction 



There are holographic dual models based on type IIA and type IIB string theories which 
describe certain aspects of strong interactions when matter fields transforming under the 
fundamental representation of the gauge group are taken into account [1, 2, 3, 4, 5]. In these 
D-brane constructions mesons are modeled as fluctuations of flavour branes in the probe 
approximation. From some of these models it is possible to derive properties of mesons 
such as ratios of masses, decay constants and the chiral Lagrangian, among others. This 
makes the alluded models interesting for investigating certain aspects of the phenomenology 
of strong interactions. Mesons derived from this kind of models are generically known as 
holographic mesons. One further important and very interesting question that one may ask 
in the context of these models is about the structure of the holographic mesons, when they 
are probed by virtual photons in deep inelastic scattering (DIS) processes. 

The DIS cross section is proportional to the Lorentz contraction of a leptonic tensor, 
entirely computed within perturbative quantum field theory, and a hadronic tensor W^ u 
which receives contributions from the strong coupling regime of QCD. The hadronic tensor 
is computed using the optical theorem, and is recast in terms of the two-point correlation 
function of electromagnetic currents inside the scattered hadron. From the most general 
Lorentz decomposition of this tensor the so-called structure functions can be extracted. 
In this work we shall be focussed on DIS from mesons. Spin-zero mesons have only two 
structure functions called Fi and F 2 . In addition, polarized spin-one hadrons have eight 
structure functions [6], as we shall review in the next section. 

From the point of view of the gauge/string duality [7, 8, 9], Polchinski and Strassler [10] 
proposed a method for deriving the hadronic tensor and the DIS structure functions for 
glueballs and spin-^ hadrons, focussing on confining gauge theories that are conformal (or 
nearly conformal) at momenta well above the confinement scale A. For instance, one can 
think of the N = 1* SYM theory [11, 12], that can be obtained from the Af = 4 SYM theory 
explicitly broken at a certain mass scale to pure M = 1 SYM theory. Further developments 
have applied the ideas of [10] to scalar mesons, unpolarized vector mesons [13, 14] and 
nucleons [15, 16]. On the other hand, DIS structure functions [17] and lepton-pair photo- 
production rates [18] have also been studied in the strong coupling regime of an M = 4 
SYM plasma using a holographic dual description to obtain correlation functions of two 
electromagnetic currents. More recently, it has been considered 0(a' 3 ) corrections from type 
IIB string theory on DIS [19], plasma photoemission [20], as well as other plasma properties 
which rely on the computation of two-point correlation functions of electromagnetic currents 
such as the electrical conductivity [21, 22]. 

In this paper we investigate DIS from polarized spin-one mesons and obtain the DIS 
structure functions by extending the proposal of [10]. We derive the interactions in the bulk 
directly from the expansion up to quadratic order in the derivatives of the fields of the Dirac- 
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Born-Infeld action of the flavour branes in the probe approximation for the D3-D7-brane 
model [1] and the D4-D8-D8-brane model [5] . The derivation we use, which is very different 
from the one presented in [13, 14], is totally general and can be applied to other holographic 
dual models and for different holographic hadrons. We also obtain the corresponding DIS 
structure functions for scalar mesons using these models. 

The main difference with [13, 14] is the way we construct the interaction in the bulk. 
So, let us summarize the method we use in this paper. The holographic dual models we 
consider have an internal Einstein manifold which contains a sphere 1 : i.e. the D3-D7 brane 
model has an S 3 and the Sakai-Sugimoto model has an S* 4 . The isometry groups are SO (A) 
and 5*0(5), respectively. We can consider a U(l) subgroup of these isometry groups. From 
the point of view of the dual gauge theory, these global bosonic symmetries correspond to 
the i?-symmetry group SU(2)ji in the D3-D7 brane model, while for the Sakai-Sugimoto 
model this corresponds to the global symmetry group SO (5). Now, consider a field in the 
probe-brane world-volume in each of these models. Its wave function can be factorized as 
the product of a plane-wave in the four- dimensional Minkowski spacetime times a function 
depending on the radial coordinate, times a spherical harmonic on the corresponding sphere. 
Notice that this spherical harmonic satisfies an eigenvalue equation, with eigenvalue equal 
to the charge under the referred global symmetry, let us call it Q. Also, recall that for each 
global continuous symmetry there is a corresponding conserved Noether's current. This is 
the situation in the bulk. Now, suppose we insert a U(l) current at the boundary theory, 
this induces a fluctuation on the boundary conditions of the bulk fields coupled to this 
current at the boundary. The fluctuation propagates in the bulk. The precise form of this 
fluctuation comes from the fact that this boundary current is a Lorentz vector which couples 
to the boundary value of a gauge field propagating in the bulk. As we shall explain later, 
this is an off-diagonal fluctuation on the bulk metric of the form 5g m j = A m (y,r)vj(Q), 
where A m is the gauge field and Vj is the Killing vector corresponding to the isometry 
subgroup U(l) C SO(n + 1) of the S n of each considered model. In this way, it is natural 
to straightforwardly construct the interaction Lagrangian in the bulk by coupling the bulk 
Noether's current above to the gauge field A m with a strength given by the charge Q. So, this 
is the main point in the construction that follows here. The existence of the bulk conserved 
Noether's current is crucial in this construction. 

The main new results we present in this paper can be summarized as follows. We develop 
a consistent method to derive the interaction Lagrangian from the Dirac-Born-Infeld action 
of flavour probe branes. Then, we obtain the full set of structure functions for scalar and 
polarized vector mesons in both holographic dual models. In addition, we find a completely 
new full tower of vector mesons from the Sakai-Sugimoto model. The reason to use these 
two very different holographic dual models with flavours to study deep inelastic scattering 

1 In a more general holographic dual model the sphere could be replaced by any other compact Einstein 
manifold, and the same argument holds. 
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from polarized spin-one hadrons is to be able to compare the structure functions and discuss 
their model-dependent as well as model-independent properties. We have obtained analytic 
expressions for the full set of eight structure functions: F±, F2, gi, g 2 , bi, b 2 , &3, &4, which 
come from the standard decomposition of the hadronic tensor of spin-one hadrons. They are 
functions of two independent dimensionless kinematical variables: x (the Bjorken parameter) 
and t, both to be defined in the next section. Indeed, we show that these functions split into 
a model-dependent factor, which is common to the eight structure functions of each model, 
and a mo del- independent one. We discuss model-independent properties of the structure 
functions, such as the relations 2Fi(x) = F 2 (x) and 2b\{x) = b2(x), which we have obtained 
for all values of the Bjorken parameter (neglecting terms proportional to t), though the 
supergravity approximation only holds for values of x close to one 2 . We have found a 
missing x-factor in comparison with the usual Callan-Gross relations. This is due to the 
fact that, within the pure supergravity description, lepton scattering is produced from the 
entire hadron [10]. Also, we have found the relation b\ (x) = 3Fx(x) in agreement with the 
expectations from [6] for the p-meson, as well as an additional set of new relations among 
some of these structure functions that we explain in the last section of the paper. 

In Section 2 we describe the kinematics of the DIS from polarized spin-one mesons which 
shall be relevant for our studies in the following sections. In Section 3 we very briefly review 
the proposal of reference [10]. In Section 4 we begin with a brief review of the D3-D7- 
brane model description following [1] and we show the expressions for the scalar and vector 
mesons. Then, we develop a consistent method to obtain the interaction Lagrangian between 
the mesons and the gauge field which arises from the fluctuation of the metric induced by 
the insertion of a current at the boundary of the spacetime. We obtain the hadronic tensor 
and the corresponding eight DIS structure functions from polarized spin-one vector mesons 
of the dual M = 2 gauge theory. In Section 5, we carry out a similar programme for 
the Sakai-Sugimoto model [5] using our method to derive the interaction Lagrangian from 
first principles. We have also obtained a full tower of vector fluctuations of the D8-brane, 
which has an additional quantum number, £, compared to the mesons considered in [5], 
that comes from the expansion on spherical harmonics on S* 4 . Besides, due to the fact that 
the DIS relevant bulk interaction region corresponds to large values of the four- dimensional 
energy, we consider the metric for that specific region. This allows us to calculate full 
analytic expressions for these mesons depending on the nine coordinates on the D8-brane. 
We extensively discuss our results in the last section of the paper. 

2 Initial gauge/string duality studies for small values of the Bjorken parameter can be seen for example 
at references [10, 23, 24, 25, 26]. 
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2 Kinematics of DIS from spin-one hadrons 



Deep inelastic scattering is a high energy process which allows to investigate the hadronic 
structure. Typically, a lepton is scattered from a hadron in a kinematical regime where the 
hadron becomes fragmented in many particles which are not measured. The process can 
be described as the electromagnetic scattering of a lepton by a quark or a parton inside 
the hadron. In this section we follow the conventions of Manohar [27], except for the four- 
dimensional Minkowski metric that we take as r]^ v = diag(-l, +1, +1, +1). So, let us briefly 
summarize the kinematics of DIS from vector and scalar mesons following the quantum field 
theory analysis from references [27] and [6]. 

Let us consider an incoming lepton beam with four-momentum (where k° = E) which 
will be scattered from a fixed hadronic target. The four- momentum of the scattered lepton 
k'^ (with k'° = E') is measured, but the final hadronic state (that we call X) is not. The 
lepton exchanges a virtual photon of four-momentum g M with the initial hadronic state. 
The virtual photon probes the hadron at distances as small as 1/ y/q 2 . If the hadron is not 
fragmented the scattering is called elastic, otherwise the resulting process is DIS. 

The DIS amplitude is given by: 

iM = He) 2 {^^j {hVm\Ksi)(x\r h {Q)\p,h), (1) 

where e is the lepton electric charge, s; is the polarization of the incoming lepton, h is the 
polarization of the hadronic initial state, and Jf, are the electromagnetic currents of the 
lepton and hadron, respectively. The polarization h can be chosen as the spin in the direction 
of an arbitrary axis, usually the direction of the incident beam. The differential DIS cross 
section is then 



da 



v ., (2tt)32£' v ' v ^>{2E){2M){v rel = l) 

v r d 3 k' (2n) 4 5 4 (k + P - k' - P x ) e 4 
Y-/ (2tt) 3 2£' (2E)(2M) ~q A 

x(P,/ i |^(0)|X)(X|jnO)|P,/i)(A;, Si |J^(0)|A;')(A;V^(0)|A;, Si ) 



(2) 



where P M and P£ are the hadronic initial and final momentum. M 2 = —P 2 and M\ = —Px 
are the initial and final hadronic squared masses, respectively. The relation (a\J fl \f3)* = 
(/3|J M |a) has been used, since = J M . Notice that since the polarizations of the final 
lepton and hadron states are not measured there is a sum over them. 
The leptonic tensor is defined as 

E (k'\J?(0)\k, Sl )(k, Sl \J?(0)\k'), (3) 

final spin 
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and can be computed pertubatively within quantum electrodynamics. For a spin-| lepton it 
is given by 

= 2[W + k v k'» - rf v {k ■ k' - mf) - ie^q a s w ], (4) 

where mi is the lepton mass. 

On the other hand, the hadronic tensor W^ v is given by [27] 

W^(P,q) h , h = i- J d 4 xe^ x (P,h'\[J^x),J u (0)]\P,h), (5) 

where h and h' are the polarizations of the initial and final hadronic states, respectively. By 
inserting a complete set of eigenstates and taking into account translational invariance, the 
above expression becomes 



W^(P, q) Wh = {2nf5\q + P- P X )(P, h'\J^0)\X)(X\J v (0)\P, h) 



4vr x 



-{2nf5\q -P + P X )(P, ti\J v (0)\X)(X\Jv(0)\P, h) 



(6) 



The allowed final states must satisfy the condition P x > P°, since M\ > M 2 . Since q° > 0, 
only the first delta function in Eq. (6) fulfills this condition, thus the sum in reduces to 
Eq.(2). Finally, one gets 

' /2 °' ' E ' -r V W» v {P,q) h , h . (7) 



dE'dtt 16vr 2 g 4 ME 

It is important to notice from Eq.(4) that the spin-independent part of the leptonic tensor 
is symmetric under //■<-)■ u, while its spin-dependent part is antisymmetric. Therefore, an 
unpolarized lepton beam only probes the symmetric part of the hadronic tensor W^ u . 

The relevant hadronic structure for DIS can be completely characterized by W^ u . The 
probability that a hadron contains a given constituent with a given fraction x of its total 
momentum is given by the partonic distribution functions, which unfortunately cannot be 
computed with perturbative QCD. This is because they depend on soft (non-perturbative) 
QCD dynamics which determines the hadronic structure as a confined state of quarks and 
gluons. Therefore, the incoming lepton gets scattered from the hadronic constituents which 
carry a fraction x of the total hadronic momentum. 

For the hadrons composed by massless partons the probability of finding a parton with a 
momentum xP^ is given by the distribution function f(x, q 2 ). When the partons are free, this 
function becomes independent of q 2 , i.e. / = f(x), which is the so-called Bjorken scaling. 
Obviously, this scaling does not hold for QCD where the partonic distribution functions 
change as a function of q 2 , since each parton tends to fragment into multiple partons with 
lower x. Thus, the hadronic structure in QCD does depend on q 2 , being the parton number 
increased while the average value of x decreases as q 2 increases. 
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Now, let us focus upon the hadronic structure functions. First, let us review a few basic 
properties. They are dimensionless functions depending on P 2 , P ■ q and q 2 . It is useful 
to write them in terms of the dimensionless variables t = A- and x = — 7^3— . The ranges 
of these variables are < x < 1 and t < 0. The structure functions are obtained from 
the most general Lorentz decomposition of the hadronic tensor, satisfying certain physical 
requirements for W^ v . For parity preserving interactions it is necessary to satisfy current 
conservation d^J^(x) = which implies that 

q'W^P, q, h) = q'W^P, q,h)=0. (8) 

It must also be parity invariant: 

W^(P, q, h) = W w (P p , qp, h p ) , (9) 

where the subindex p indicates parity-transformed quantities. Time-reversal symmetry must 
also be satisfied, which means that 

W^(P, q, h) = W: TtXT {P T , q T , h T ) , (10) 

where the subindex T means time reversal. In addition, there is an identity relating for 
q° > with W^ u for q° < 0: 

W lu ,(P,q,h) = -W vli (P,-q,h), (11) 

which comes from translation invariance. At this point we introduce the hadronic tensor for 
polarized spin-one hadrons for parity-preserving interactions as derived in [6] by Hoodbhoy, 
Jaffe and Manohar. There are four new structure functions in comparison with spin-| 
hadrons, which are called bi, b 2 , b^ and 6 4 . They contribute to the symmetric part of W^ u , 
thus contributing to DIS from an unpolarized target. Therefore, the hadronic tensor can 
be written in terms of eight independent structure functions. Omitting terms containing 
and q u , since they vanish after the contraction with the leptonic tensor, the hadronic tensor 
reads 

F2 b 2 &3 

W^ u = F x r]^ v - — -P^Pu + &i?> - — {s^ + V + - ~2 (^V - *v) 

&4, , x igi \ a W2 A/ D a ™\ 

~^{s^ - V) - -p7^e Mi , Aa g s - - ^ 2 e^\aq {P • q s - s-qP ) , 

(12) 

where we have used the following definitions: 

= (p l g)2 (g • C* g • C - ^P • g) 2 «)^ , (13) 
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V = (p^p ( 9 • C * q • C - V P • 9) 2 ") P /* P " , (14) 

(15) 

*V = (QCv + + - \p,Pu) , (16) 



-.p^CC^p, (17) 

being k = l—Ax 2 t and s a a four- vector analogous to the spin four- vector for spin-| particles. 
Besides, ( and (* denote the initial and final hadronic polarization vectors, respectively. The 
condition P ■ ( = is satisfied, and the normalization is given by ( 2 = —M 2 . 

Recall that DIS amplitudes can be obtained from the imaginary part of the forward Comp- 
ton scattering amplitudes. In particular, from the matrix element of two electromagnetic 
currents inside the hadron, the tensor T^ u is defined as 

= i{P, Q\f(J,(q) J„(0))|P, Q) , (18) 

where, as before, P is the four-momentum of the initial hadronic state (where for brevity 
we have omitted its Lorentz index), q is the four- momentum of the virtual photon and Q 
is the charge of the hadron. T(Oi0 2 ) indicates time-ordered product between 0\ and 2 
operators. The tilde stands for Fourier transform. The tensor = T M „(P, q, h) has identical 
symmetry properties as W^ u (P,q,h), thus having similar Lorentz-tensor structure to W^. 
Using the optical theorem one gets 

ImFj = 27rFj, (19) 
being Fj the j-th structure function of T^ u and Fj the corresponding one of W^ v . 



3 A holographic dual description of DIS 

Let us very briefly describe the idea presented in [10] by Polchinski and Strassler to study 
DIS for gauge theories that have holographic dual descriptions. Full details are explained 
in the original reference. Within the supergravity approximation they calculate hadronic 
structure functions for x ~ 1 (with x < l) 3 . In particular, they consider confining gauge 
theories in four dimensions such as certain deformations of Af = 4 SYM, from which they 
study DIS from glueballs and spin-| hadrons. The examples considered in that paper are 

3 In [10] also a small- a; calculation was carried out in terms of a string theory analysis, i.e. beyond the 
pure supergravity approximation. 
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UV conformal or nearly conformal. Thus, the dual string theory is defined on the AdS 5 x W 
background, where W can be an Einstein manifold. The metric is given by 



r 2 R 2 —2 

ds 2 = -^V^dy^dy" + —dr 2 + R 2 ds w , (20) 

where R = (ATig s N) l ^a a / 2 is the AdS 5 radius when W is the five-sphere 4 . The four coordi- 
nates y^ are identified with those of the gauge theory, while r is the holographic radial one. 
Coordinates on W are denoted by Q. The ten-dimensional energy scale is given by (up 
to powers of the 't Hooft coupling A = gy M iV = 47rg s N, with the string coupling g s ), while 
the characteristic four- dimensional energy is 

£ (4) ~^- (21) 

If one is interested in the large N limit of confining gauge theories, the geometry of the 
holographic dual model at large r is approximately that of Eq.(20), however, it must be 
modified at a radius corresponding to 

r ~ r = AR 2 , (22) 

where A is the confinement scale. Nevertheless, the dynamics of interest for q ^> A corre- 
sponds to r int ~ qR 2 ^> r , where the conformal metric (20) can be used. r int denotes the 
bulk region where the relevant interaction occurs as we explain below. 

Polchinski and Strassler employ the dual string theory description to obtain the matrix 
element T^ v . Its imaginary part is given by 

ImT^ = 7T ^2(P,Q\r(0)\P x ,X)(P x ,X\J»(q)\P,Q), 

Px,X 

= 2tt 2 £ 5(M 2 X + [P + q] 2 )(P, Q\r(0)\P + q, X)(P + q, X\J"(0)\P, Q) . (23) 

x 

The momenta P M , g M , the polarization n M and the currents J M are considered as four- 
dimensional quantities, being their Lorentz indices raised and contracted with rf v . 

In the large N limit of the gauge theory only single hadron states will contribute. In the 
gauge theory for the s-channel we have 

S = -(P + g) 2 ^g 2 Q-l), (24) 

where it has been used —P 2 <^ q 2 . The corresponding ten-dimensional scale is 

S = - 9 -p„,p„ < -r(p + ,wp + - £ e (i - 1) = ■ (25) 



4 Otherwise it must be other volume factors multiplying the relevant factor (<7 s 7V) 1/4 o:' 1/2 of the metric 
scale, however, for those cases the present discussion also applies. 



8 



The 't Hooft parameter appears in the denominator, so if {g s NY l l 2 « i< 1 we have 
a's <C 1. Therefore, for x ~ 1 (with x < 1) only massless string states are produced, and we 
are dealing with a supergravity process [10]. This is the regime studied in the present paper. 

Now, let us focus on the way that the DIS process is viewed from the bulk theory. Recall 
that in the four-dimensional boundary theory we have to obtain the two-point function of 
two electromagnetic currents inside the hadron. Using the gauge/string duality the current 
operator inserted at the boundary of the AdS space induces a perturbation on the boundary 
condition of a bulk gauge field. The perturbation excites a non-normalizable mode which 
propagates within the bulk [8, 9]. We can see how bulk gauge fields appear as follows. 
The isometry group of the manifold W corresponds to an R-symmetry group on the field 
theory. We can consider a U(1)r subgroup and the associated R-symmetry current can now 
be identified with the electromagnetic current inside the hadron. In addition, for the global 
symmetry group corresponding to the isometry of manifold W, there is a Killing vector Vj. 
It excites a non-normalizable mode of a Kaluza-Klein gauge field, 

5g mj = A m (y, r)vj(Q) . (26) 

In the case of glueballs, the holographic dual field in [10] corresponds to the dilaton. Thus, the 
incoming bulk dilaton field $j couples to the bulk {/(l)-gauge field (induced by a current 
inserted at the boundary) and to another dilaton <&x, which represents an intermediate 
hadronic state (See Figure 1). The latter propagates in the bulk and couples to an outgoing 
dilaton $j (corresponding to the final hadronic state) and a gauge field A u in the bulk 
which comes from the insertion of a second boundary theory current. This Witten diagram 
represents the holographic dual version of the optical theorem [10] for a DIS process which is 
schematically shown at the boundary slice, where the leptonic currents are indicated with /, 
and they couple to the virtual photons indicated by dashed lines at the boundary slice. Then, 
these photons couple to the currents inside the hadrons, which are indicated by solid lines 
on the boundary slice and are labeled by the momenta carried by the initial, intermediate 
and final hadronic states, p, Px and P/, respectively. The present holographic approach 
corresponds to a DIS process where the lepton is scattered from an entire hadron, which 
becomes excited but does not fragment. In the next two sections we extend these ideas to 
other holographic dual models. 

4 DIS from mesons in the D3-D7-brane model 

In this section we firstly give a brief introduction to the D3-D7-brane model developed in 
reference [1]. In [28] it was shown that by introducing Nf D7-branes in the AdS§ x S 5 
background of type IIB string theory, it is possible to describe Nf hypermultiplets in the 
gauge theory preserving M = 2 supersymmetries in four dimensions. The starting point is 
a set of coincident parallel N D3-branes and Nf D7-branes which share directions 0-3 with 
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Figure 1: Schematic representation of the AdS/CFT prescription for the calculation of the tensor 
T^ u . The interaction in the bulk occurs at ri nter action, we also show the slices corresponding to the 
IR cutoff tq and the boundary. 
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the former ones. The hypermultiplets in the gauge theory arise from the lightest modes of 
the fundamental strings extended between D3 and D7-branes: those are modes of the type 
3-7 and 7-3, whose mass is m q = L/2na', where L is the distance between D3 and D7-branes 
in the (8, 9) plane. In the decoupling limit of the D3-branes, namely g s N ^> 1 where g s 
is the string coupling, the background becomes AdS$ x S* 5 . Also, if N ^> Nf then the 
backreaction of the D7-branes can be neglected. In these limits the description corresponds 
to Nf probe D7-branes in the referred AdS 5 x S 5 geometry. We consider only the Abelian 
case (Nf = 1), which implies the existence of a global U(l) symmetry group, associated 
with the flavour symmetry of the M = 2 SU(N) SYM theory. This theory has dynamical 
quarks, thus it is very interesting to investigate the DIS structure of scalar and vector mesons 
which arise from this model. We should be aware that this model does not allow to describe 
spontaneous breaking of chiral symmetry. This makes a difference with respect to QCD and 
the Sakai-Sugimoto model [5]. In [1] the mass spectra for different types of mesons have 
been computed through its holographic dual description. In this work we shall obtain the 
structure functions of the hadronic tensor for DIS from those scalar and polarized vector 
mesons. 

We develop a systematic method to derive the interaction Lagrangian C int in the holo- 
graphic dual theory, in a complementary way with respect to the proposal of Polchinski and 
Strassler [10]. We emphasize that the interaction Lagrangian obtained using both methods 
is the same for the case studied by Polchinski and Strassler [10], while it gives very different 
results for vector mesons in comparison with [13] and [14]. The present approach guarantees 
the existence of a conserved current in eight dimensions straightforwardly derived from the 
Dirac-Born-Infeld action. From this current we shall see that the property q^W^ = is 
satisfied. Moreover, we have checked that Ci n t is invariant under parity and time-reversal 
transformations. Thus, the hadronic tensor automatically satisfies the properties indicated 
in Eqs.(8)-(11). All these are consistency checks that confirm the correctness of the proce- 
dure we develop here. Next step will be the derivation of the structure functions for scalar 
and polarized vector mesons, obtaining analytical expressions for all of them. 



The D3-D7-brane model 

For this brief subsection we follow reference [1]. Let us consider a type IIB string theory 
background consisting of N coincident D3-branes. This corresponds to the geometry AdS$ x 
S 5 whose metric is 

ds * = ^ ds 2 (E^) + %dZ ■ dZ , (27) 
R 2 r 2 

where the coordinates Z l , i = 1,...,6 parameterize the transversal directions to the D3- 
branes and r = \Z\. In addition, the AdS 5 radius is given by R 2 = ^4irg s Na'. 
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Let us add a probe D7-brane at a certain distance \Z\ — L from the D3-branes in the (8, 
9) plane. In this case the hypermultiplet becomes massive, and the R-symmetry is SU(2)r. 
The induced metric on the D7-brane is 

^ = l^ iAE ™ ) + ^ dp > + J^ d nl (28) 

where p 2 = r 2 — L 2 and fl 3 are the spherical coordinates in the space spanned by the 
coordinates 4-7. For L = the metric reduces to AdS 5 x S 3 , rendering a conformal gauge 
theory. On the other hand, if L ^ 0, the above metric is only asymptotically AdS 5 x 
S* 3 for p ^> L. This shows the explicit breaking of conformal invariance induced by the 
hypermultiplet mass m q = L/2na', which is recovered in the high energy limit, E ^> m q . 
Notice that the radius of the S 3 is not a constant, and moreover it goes to zero for p = 
(which corresponds to r = L), where the D7-brane ends as seen from the projection on the 
AdS 5 . 

In [1] the spectra of scalar and vector mesons were computed, and they were arranged 
in supersymmetric multiplets which transform according to representations of the global 
SU{2) R x SU{2) L symmetry group. These mesons correspond to excitations of open strings 
of the D7-brane. The dynamics of the probe D7-brane fluctuations is described in terms of 
the following action 

Sdi = -P7 j d 8 ^-det(P[g] ab + 27ra'F ab ) + ^^p 7 J P[C {4) ] AFAF, (29) 

g ab stands for the metric (28), p 7 = [(2k)' 1 '^q;' 4 ] -1 is the D7-brane tension and P denotes 
the pullback of the background fields on the D7-brane world-volume. The relevant part of 
the Ramond-Ramond potential in the Wess-Zumino term is 

4 

C (4) = ^dy° A dy 1 A dy 2 A dy 3 . (30) 

Next, in order to fix our notation and give a self-contained derivation of the hadronic 
tensor, we show briefly how the mesons are obtained from the present model. Then, we 
carry out the derivation of the hadronic tensor and the structure functions for scalar and 
vector mesons. 

DIS from scalar mesons 

The equations of motion for scalar mesons are obtained from transversal fluctuations of the 
D7-brane [1] 

Z 5 = + 2na'x, Z % = L + 2ira'ip , (31) 
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where the coordinates Z 5 and Z 6 lie on the (8, 9) plane, transversal to the D7-brane. x an d 
if are the scalar fluctuations whose Lagrangian is straightforwardly derived from action (29), 



C = -p 7 ^\det(g ab )\^det(5 ab + (g ac R 2 /r 2 )(2rfa') 2 (d cX d b x + d c tpd b tp)) . (32) 

All indices denote directions in the world-volume of the D7-brane. Expanding the above 
Lagrangian up to quadratic order in the derivatives of the fluctuations, the previous equation 
becomes 



C = -p 7 J\det(g 



Jab 



2R 2 

1 + -^r9 a \™f{d bX daX + dwd a <p) 



(33) 



which only depends upon derivatives of the scalars. The D7-brane wraps the S* 3 , whose radius 
is denoted by p. Substituting r 2 = p 2 + L 2 and the metric (28) in the quadratic Lagrangian 
one obtains the equations of motion (EOM) for scalar fluctuations of the D7-brane in the 
probe approximation, which reads 



*($^**)= ' (34) 

where <£> is either fluctuation (x or ip), while is the metric of a 3-sphere of unit radius, 
which together with p include the directions (Z 1 , ■ ■ •, Z 4 ). The EOM can be written more 
explicitly as 

L> /; -S^d^Q + d p ( - 2 P \ 2 , d ( A + / 3 T2 y Vi<t> = , (35) 



(p 2 + L 2 ) v P \(p 2 + L 2 ) p ) (p 2 + L 2 ) 

where V« is the covariant derivative on S* 3 . Indices a, b, c, • • • denote coordinates on the 
D7-brane world-volume; i, j, k, ■ ■ ■ are on S 3 ; p, u, ■ ■ ■ indicate directions parallel to the 
D3-branes and m, n, ■ ■ ■ correspond to 0, 1, 2, 3, p coordinates. 

Now, we need the wave functions for the scalar mesons in a region where 



/■,■„, = VPL + L 2 ~ qR 2 > r = AR 2 > V2L p int > L , (36) 

where rj nt denotes the interaction region and A is an infrared cutoff for the four- dimensional 
gauge theory. In this region the metric can be approximated by AdS$ x S* 3 , which reflects 
the fact that for energies as high as E ^> m q the explicitly broken conformal invariance gets 
restored. Thus, the resulting EOM is 

y (P'djAS + d p (p%$) + pV'V^ = . (37) 

The proposed Ansatz for the solution is 

& = (f)(p) e ip - y Y e (S 3 ) , (38) 
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where Y E (S 3 ) are the scalar spherical harmonics on S* 3 , which satisfy the eigenvalue equation 

V i V l Y i {S 3 ) = -1(1 + 2)Y\S 3 ) . (39) 

Notice that in order to simplify the notation, from now on we will call $ = $ £ , omitting the 
£ index. 

Thus, the Lagrangian (33) for the complex scalar <3> in the interaction region can be written 

as 

C = -p 7 (na') 2 ^tg~\^g ab (d a ^d b ^ + d b ^d a ^) . (40) 

V pZ 

Now, we have to derive the EOM for the fields in the AdS space. Thus, replacing the solution 
(38) for the scalar mesons in the interaction region defined by Eq.(36), we obtain 

d 2 p <j) + ^ l d p (t) + (-p 2 /?V 4 - ™?P~ 2 )0 = o , (41) 

with m\ = £{£ + 2). 

For the initial/final hadronic state (IN/OUT), with four- momentum P, we can use the 
leading behaviour of the solution in the region p ~ p int , while for the intermediate state X 
with four- momentum Px, we have to use the full solution. Thus, we obtain 

^in/out = Q(Ap)- l - 1 e^y(Q) , (42) 
$ x = Cx K-VH^p-\j u {^^y p - yY{n) , (43) 

where J u is the Bessel function of first kind, v = \Jm 2 + 1 = £+1 and s = —(P + q) 2 = M x is 
the mass-squared of the intermediate state, while cx and q are dimensionless constants. In 
order to study the parity we need to use the property [F £ (f2)]p = (—1) Y (JT). As a result, 
the mesons are scalar and pseudoscalar mesons for even and odd values of £, respectively. 

At this point we have all the ingredients to calculate the hadronic tensor using the holo- 
graphic dual prescription. For this purpose we first consider that the holographic meson 
couples to a gauge field in the bulk of the string theory dual model. In order to understand 
this coupling let us consider the insertion of a £/(l)-current at the AdS 5 boundary. This 
induces a perturbation on the boundary conditions and excites a non-normalizable mode A m 
in the AdS^, interior. This mode corresponds to the extension towards the AciS'-interior of 
the A^l^d field, which is in fact its boundary value. The perturbation propagates inside the 
AdS, inducing a metric fluctuation which is the product of the gauge field A m (y,p) and a 
Killing vector Vj corresponding to an isometry of the internal space S* 3 . Thus, we consider 
the fluctuation proposed in [10]. Perturbing the metric in Eq.(40) we have g ab g ab - 5g ab , 
where 5g ab = A m (y , p)v' J \Q) and using v j djY e (n) = iQ e Y e (Q) (below we shall drop the £ 
label from the spherical harmonics and from Qg), the interaction Lagrangian we obtain is: 

C int = tQp 1 ^a') 2 ^ig~\R 2 p^ 2 A m ^d m ^ x - <^<9 m $) . (44) 
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The angular dependence on the spherical harmonics corresponds to functions which are 
charge eigentates, with charge Q under the U(l) symmetry group induced by transformations 
on the internal S* 3 in the direction of the Killing vector vK 

Alternatively, Ci nt can be obtained from the coupling of the gauge field A rn to the Noether 
current corresponding to the internal symmetry of the action (40) under global phase trans- 
formations associated with the U(l) group, which is a subgroup of the 5*0(4) isometry group 
of S 3 , generated by a Killing vector v j . The Noether's current is: 

r = aJ5f) 6 + m^) e ' ■ (45) 

Now, for an infinitesimal transformation with parameter i?/2(= constant) we have 

5$ = — «-$ =4> = — -$ , 

S$* = it§* e* = l$*. (46) 
2 2 v ; 

Therefore, the Noether's current reads 

j m = i(x 7 (na') 2 R 2 p- 2 (<$>d m <S>* x - $^<9 m $) . (47) 



If we define Ci nt = Q^J\detg\A m j m , we obtain the same Ci nt given by Eq.(44) from a metric 
fluctuation. Consequently, Ci nt is given by the coupling of the gauge field A m to the conserved 
Noether's current j m . Now, a crucial point for the present derivation of the hadronic tensor 
is that the fields (38) are charged under a U(l) symmetry group, which is a subgroup of the 
5*0(4) isometry group of S 3 , generated by a Killing vector vK In order to check this we may 
ask how they transform under a coordinate transformation corresponding to an isometry 
on S* 3 : (x m , x l ) — > (x m , x % + v l e(x m )), where v % is the Killing vector which generates the 
isometry. We consider i > which means that the scalar fields are charged under a U(l) 
group. 

In order to calculate the matrix element relevant for the hadronic tensor, we will use the 
identification [10] 

S int = (2tt)V(Px -P- q)n,(P + q, X\ J"(0)|P, Q) . (48) 

Therefore, we must calculate Si nt . We start by deriving the coupling A m j m , using the 
expressions 

P 

A p = -- 2 ^q,d p A U) (49) 



A, = n^KJ^), 
P V p J 
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where and A p are the components of A m (y,p) deduced in [10]. The current conservation 
equation is d u j u + p~ 3 d p (p 3 j p ) = 0. We obtain 

A m j m = A^f - i^a^ 7 ) - i^d p {A^f) . (50) 

Then, the interaction is given by 

S int = / d 8 xQ^tg\A m j m 

Jro 

= d*x2p 7 (irafQj\a^\p- l R*(p» + £)A^*x 

/ A ^Qy1det^(^p 3 J 0|^ • (51) 

The integral can be calculated with the substitution z = i? 2 /p and taking 1/A — > oo since 
goes to zero very fast for p < pi nt . We obtain that the second term of Eq.(51) evaluated at 
the integration limits vanishes. After integration of the remaining term, and using Eq.(48), 
it results 

{P + q,X\J»(0)\P,Q) =« 1 / 1 (p/' + £), (52) 

where K\ = 2pjQ{jia') 2 CiC* x A u ~ 5 ^ 2 q 3 ^ 2 x~ 1 ^ 4 (l — x) 1 ^ and the result of the integral in p is I\. 
For \t\ < 1 we can approximate s ~ q 2 (l/x — 1). Doing this, we obtain 

h = 2 1+ T(2 + u)q-^ u (l - x) u/2 x u/2+2 . (53) 

In order to compute ImT M/3 as in [10], we have to multiply Eq.(52) by its complex conjugate 
and sum over radial excitations. For this we have to know the density of states, which can 
be estimated by introducing an IR cutoff at r . In this way the distance between zeros of 
the Bessel function of Eq.(43) gives M n = nnA. In the large N limit this becomes a sum of 
delta functions, and for large q we have 

£«-*)- (^) _1 ~ (2,^-1. (54) 

Assembling all these partial results, we obtain 

ImT^ = 47r 5 Q 2 2 2+ ^[r(2 + z/)]V 7 («') 4 |c J | 2 |c x | 2 

xA 2 "-V 4 - 2 "(l - xyx^(p» + £) (V + . (55) 

The remaining step is to obtain the structure functions for the scalar mesons from Eq.(19). It 
can be seen that satisfies the properties described above. Then, the structure functions 
we obtain are: 

/// 2 r/ 4 \ /\ 2 \ u+l 

Fi = o, F2 = A ° Q i^){^) ^ +3 a-*r> ( 56 ) 

where A' = 47r 5 2 2+2i/ |ci| 2 |cx| 2 [r(2 + u)] 2 is a dimensionless normalization constant. We shall 
comment about these results in the last section of the paper. 
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DIS from vector mesons 

Now, we focus on how to derive the structure functions for polarized vector mesons in the D3- 
D7-brane model. First we derive C int for the coupling of a gauge field A m to vector mesons. 
For this purpose we must find explicit expressions for the fields. With this we obtain all 
the structure functions for polarized vector mesons, which have not been considered in the 
previous literature in the context of the gauge/string duality. 

In [1] the spectrum of vector mesons has been computed. They arise from fluctuations 
of the vector fields of the Dirac-Born-Infeld (DBI) action of the probe D7-brane, in the 
directions parallel to this brane. The starting point is the action (29). As before, first the 
EOM are derived from it and then the Lagrangian is expanded up to quadratic order. This 
new Lagrangian gives the same EOM as (29). Recall that in (29) F ab = d a B b - d b B a . The 
EOM is 

d a (^tg~\F ab ) - ^p(p 2 + L 2 )e bjk djB k = , (57) 

where e^ k is the Levi-Civita pseudo-tensor, the indices a, b, c, d, ■ ■ ■ run over all directions 
of the D7-brane world-volume, and i, j, k, ■ ■ ■ belong to S 3 . Also, the second term, coming 
from the Wess-Zumino Lagrangian, is only present if b is on S* 3 . The Ansatz [1] for the 
solution of the vector mesons is: 

B^C^(p)e iP - y Y e (S 3 ) : P-( = 0, B p = 0, B, t = , (58) 

where there has been done an expansion in spherical harmonics on S 3 , <p(p) is a function to 
be determined, C M is the polarization vector and the relation ( ■ P = comes from d^B^ = 0. 

The interaction region corresponds to pi nt ^> L with L / 0, where the energy in four 
dimensions specified in (21) is A. As in the case of the scalar mesons above, in this 

regime the metric can be approximated by a pure AdS$ x S* 3 metric. Now, in the interaction 
region of the AdS space the EOM for b — p: is given by 

1 d m ^\detg AdS \d m B») - ^B» = , (59) 
yJ\detg A ds\ K 

where gAds is the AdS$ metric, m, n, ■ ■ •, denote 0123p coordinates and mj = £(£ + 2). 
Replacing the proposed Ansatz (58) in Eq.(59) gives the same equation as Eq. (41), with 
the solution 

c /s 1/2 R 2 \ 

4(p) = l Jv (?—?Ly (60) 

where J u is the Bessel function of first kind, c is a dimensionless constant and v — £ + 1. For 
the initial and final hadronic states we can consider only the leading behaviour for \P 2 \ <ti q 2 , 
so the solution is 

B, in/ out = C„c i A- 1 (Ap)-"- 1 e <p -» Y\S 3 ) , (61) 
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where q is a dimensionless normalization constant. For the intermediate state we consider 
the full solution 

-1/4 A -3/2 ,1/2 02. 

= Cx,c x S A p J v (LJLyiPx-v Y \S 3 ) , (62) 

where cx is a dimensionless normalization constant. In this case the radius of the Bessel 
function s l l 2 R 2 is much larger than r = AR 2 , implying that the full solution must be 
used. To study the parity transformations, using the relations [F £ (f2)]p = (— l) £ F £ (f2) and 
[C/Jp — (C°> — C) we can classify the mesons as vector mesons for even values of £ and axial 
vector mesons for odd values of i. 

From the expansion of in spherical harmonics on S* 3 , it can be seen that the gauge 
fields on the D7-brane correspond to charged fields in AdS$. To see this we can write the 
full Kaluza-Klein expansion 

B, = H^Bl (63) 

e 

where at are coefficients. Then, we consider coordinate transformations on the higher- 
dimensional spacetime such as (x m ,x l ) — > (x m ,x l + v l e(x m )), where v % is a Killing vector on 
S 3 , i denotes S* 3 coordinates and m belongs to AdS$. For the variation of the vector field 
under this transformation (being a gauge field on AdS$ x S* 3 ) we find 

SB^ = -eix^v'diBp - d^{x m ) , (64) 

using the eigenvalue equation djY £ (S 3 ) = iQeY e (S 3 ), 

SB° = -<9 M e(x m ) ^ = 0, (65) 

5Bl = - e (x m )iQ e B^ £>0. (66) 

Thus, we obtain a gauge field B® since it transforms accordingly under the U(l) transforma- 
tions parameterized by e(x). On the other hand, all the fields with £ > have a charge Qt, 
since they transform with a phase under the U(l) induced by a transformation on S 3 in the 
direction of the Killing vector. In addition, all the B e with £ > are massive vector fields 
from the AdS 5 point of view, being their masses given by m 2 = £(£ + 2) (the dimensional R 2 
factor has been ignored), since they are not invariant under gauge transformations. Their 
EOM is Eq.(59) 5 . The EOM for the vector mesons in the interaction region can also be 
derived from the following quadratic Lagrangian, 



2MW)V|det£|F ab F a V (67) 



5 We will use = as our notation, with I > 0, therefore there is a field B^ for each i. 



18 



Now, let us focus on the bulk interactions. Considering the metric fluctuation 5g m j = 
A m (y,r)vj(Vt) and v j djY(£l) = iQY(Vt) we obtain the following interaction Lagrangian, 



dm = 2tQfi 7 (7ra'yJ\detg\A m [B* Xn F nm - B n (F 



(68) 



which alternatively can be obtained from the coupling of the gauge field A m to the Noether's 
current corresponding to the internal symmetry of the action (67). This corresponds to U(l) 



transformations, thus £j nt = Qy \detg\A m j m where 

j m = 2zM^') 2 [B* Xn F nm - B n (F 



(69) 



Therefore, Ci n t results from the coupling of the gauge field A m to the conserved current j m . 
At this point, we need to compute the coupling A m j m . We use a similar method to that 
described before for scalar mesons, thus 



Sint = [°° d 8 xQj\detg\A m f 

Jrn 



TO 
COO 



= £a2 M7 K) 2 q v S^VV' 1 ^-^ I (70) 



where = (^B and 



N» = 2(C • Cx) (p" + £) + (Cx ■ q)C - (C " q)Cx ■ 



(71) 



We consider z = R 2 /p and 1/A — > oo in order to calculate the integral analytically, for 
the same reason as in the previous case. We obtain again that the second term of Eq.(70) 
evaluated at the integration limits vanishes. 

Then, comparing with the Ansatz in Eq.(48), we find 

(P + q,X\r(0)\P,Q) =k 2 I 2 N^ (72) 

where k 2 = 2p 7 Q(-na') 2 CiC* x q 1 l 2 h v ~ 7 l 2 x 1 l i (l — rr)^ 1 / 4 and the integral in p is J 2 , using s ~ 
q 2 (l/x — 1) this reduces to 

I 2 = 2^ +1 i> + 2) g—V/ 2+2 (l - xf/ 2 . (73) 

In order to obtain lm.T tJ,u from Eq.(23), we have to multiply Eq.(72) by its complex conjugate 
and sum over the radial excitations and over the polarizations of the final hadronic states 
( x . The density of states can be estimated as for the scalar mesons. Then, the resulting 
expression for the tensor is 



1TK 2 K* 2 ll 

As 1 / 2 
txk 2 k\1\ 



A 



As 1 / 2 



E 



2(C • C 



x, 



x 



2(C*'Cx) P" + 



2x 



+ (C x ■ q)C - (C* • q)Cx 



(74) 
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For the solution (58) of vector mesons, we have P • ( = 0. Thus, there are three allowed 
polarizations: A = 1, 2, 3. The normalization for the polarizations is (^(Px, A) • Q(Px, A') = 
-M 2 x 5 x ,y then 



J2(x„(Px,\)Cxv(Px,\) = M xV , u + P x ,P Xu 

A 

= M 2 xVlIU + (P„ + g M ) (P v + q v ) . 



(75) 



We must sum over the polarizations of the final hadronic state X. lmT^ u can be written in 
terms of a symmetric and an antisymmetric part, also we can neglect terms proportional to 
g M and q u , thus obtaining 

]T ImTV = H, u = H% + H% , (76) 

A 

where H^ u and H£ v are the symmetric and antisymmetric parts of H^ u respectively. They 



result 



[XV 



H b „, = 



As 1 / 2 



x 1 



AP 2 {P + qf + {q.Q{q.e) 



'P-qf-P 2 q 2 



+ (P»C + P»Q(( ■ Qh 



P-q + q' 



+ (P^ + P u Q(C-q)- 



P-q + q 2 



,(77) 



and 



2nK 2 K* 2 i 2 r i 



AsV2 



2 (CuC- CuCn) 



(P ■ q) 2 - P 2 q 2 



HP,Q-P,C){C-q)- 2 

+(P,(v - PuC,)(C -q)\ 



4 p2 + 7pq + 3q 2 



AP 2 + 7P.q + 3q 2 



(78) 



Now, we have to rewrite the hadronic tensor for spin-one hadrons from Eq.(12) in order 
to compare it with H^ u . 

Using Eq.(19) and comparing it with the hadronic tensor for polarized targets of spin one, 
we obtain a system of 6 equations and 8 functions to determine, which are the structure 
functions. However, since these functions do not depend on the polarizations, there are two 
additional equations, and therefore the structure functions become completely determined. 
The structure functions can be written in terms of dimensionless variables t — < 



and x 



, as defined in Section 2, and also in terms of the model dependent set of 



2P-q 



parameters such as A, \cx\ and Q. In this way, we obtain the following set of DIS 
structure functions: 



Fi 



A(x) 



1 

12x3 



[l-x-2xt- Ax 2 t + AxH + 8x 3 t 2 ) , 



(79) 
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F 2 


= A(x 




= A(x 


b 2 


= A(x 


h 


= A(x 


h 


= A(x 


9i 


= A(x 


92 


= A(x 



where A(x) is given by 



-^-(1 - x + I2xt - UxH - 12x 2 t 2 ) , 
ox 6 



Ax 3 
1 

2X 3 
1 



24x 3 
1 



12x 3 
1 

8x 4 



[1 — x — xt) , 
[1 — x — x 2 t) , 
(1 - Ax + 8x 2 t) , 
(-l + 4x - 2x 2 t) , 



(80) 
(81) 
(82) 
(83) 



16x 4 



-3 + 3a; + Axt + 5x 2 t - 6xH - 8xH 2 ) , 
(3 - 3x - Axt + 2x 2 t) . 



A(x) = A Q 2 



A 2 

— ) x» + \l-x) 



v-l 



A 8 / \ q< 

and A = \ci\ 2 \c x \ 2 2 4+2u [T(2 + v)] 2 ir 5 is a dimensionless normalization constant. 



(85) 
(86) 

(87) 



5 DIS from mesons in the D4-D8-D8-brane model 

In this section we extend the ideas of the previous one to the model presented in [2]. This 
model is based on a D-brane construction in type IIA string theory consisting of a large 
number of D4-branes and Nf D8-D8-brane pairs in the probe approximation (though we will 
focus on the Abelian flavour symmetry only). Direction 4 of the D4-branes is compactified on 
an S 1 , where antiperiodic boundary conditions are imposed for fermions in the gauge theory, 
thus breaking all the supersymmetries. This model geometrically realizes the U(Nf) L x 
U(Nf)n chiral symmetry. The radial coordinate U, which is transversal to the D4-branes 
has a lower bound Ukk- The radius of S* 1 shrinks to zero as U — > Ukk- The D8 and the 
D8 branes merge at U — Uq > Ukk- in the resulting D8-brane there is only a U(Nf) factor 
left. This is understood as a holographic representation of the spontaneous breaking of the 
U(Nf) L x U(Nf) R chiral symmetry down to U(Nf) v . 



The D4-D8-D8-brane model 

We begin with the metric for the described system, details are given in [5], 

ds2 =Q M»dy» + f(U)dT 2 ) + Q (^- )+ U 2 dtl 2 ), 

2nN r 



'[A 3/4 

9s[ .r) 



F 4 = dC, = 



V A 



: Q , f(U) = 1 



U KK 

u 3 



(88) 
(89) 
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where y^ (/i— 0,1,2,3) and r are the D4-brane directions. The line element on the four- 
sphere is dVt\, the volume form is 64 while the S* 4 volume is V4 = 87r 2 /3. Parameters R and 
Ukk are constants. R is related to the string coupling g s and the fundamental length l s by 
R 3 = ng s Nl 3 . 

The ^-coordinate has dimensions of length and is the radial coordinate in the 56789 
directions which are transverse to the D4-branes, U > Ukk- 

The induced metric on a D8-brane embedded in the above background is 

/77\3/2 I-/7TN3/2 //?\ 3 / 2 TJ' 2 1 //?\ 3 /2 

ds 2 =Q V^df+[Q W) + (§) JJ^]dr 2 +(§) U 2 ^ (90) 
with U' = jpU. The D8-brane action is proportional to 



S m oc J d 4 ydre 4 e-^\det(g m )\ oc J d A ydrU 4 ^ f(U) + (J ^ ■ (91) 

Since the integrand does not explicitly depend on r it is possible to express r as a function 
of U and thus the metric (90) becomes 

/TTx 3/2 //?\ 3 / 2 U 8 //?\ 3 / 2 

Now, let us recall that we want to study DIS from scalar and polarized vector mesons in 
the context proposed in [10]. We shall obtain the DIS structure functions for x ~ 1 (with 
x < 1). The fundamental characteristic of the geometry (92) is the gravitational redshift, i.e. 
the warp factor multiplying the flat four-dimensional metric. The momentum = — i<9 M is 
identified with the one in the gauge theory. As seen by an inertial observer the momentum 
on the D-brane is 

1 . . 

Pn = —n=i p » ■ ( 93 ) 
\9tt\ 



The characteristic energy scale on the D-brane is R 1 , thus the four-dimensional energy scale 
is 

773/4 

E (i) ~ ^tj • (94) 

The dynamics of interest for DIS corresponds to the limit g> A, where A is the confinement 
energy scale of the gauge theory. Thus, we shall consider the interaction in the ultraviolet 
limit, being the interaction region given by 

U int ~ q 2 R 3 > Uq = A 2 R 3 = U KK ■ (95) 

In this region the metric (92) can be approximated as 

-/Jx3/2 /i?\ 3 / 2 



ds 2 



|) V^'W'+lJ) dU* + R^U^dQl . (96) 
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We shall study the matrix elements of two electromagnetic currents in the hadron in order to 
extract the structure functions using the optical theorem. Also, it is important to recall that 
as in the case of D3-D7 model, the current inserted at the boundary induces the excitation 
of a non-normalizable mode propagating in the bulk. The boundary conditions for A m are 
given by 

KmA^U) = n^y = A,(y)U d , (97) 
lim Au(y,U) = 0. (98) 

U— >oo 

The gauge field satisfies the Maxwell equation D m F mn = 0, in the space given by 0123C7- 
coordinates with the metric (96). It is convenient to use the following gauge: 

R- 3 U 3/4 du(U 9/4 Au) + iq v A v = . (99) 

We propose the Ansatze 

A, ~ n^mU) , (100) 

A v ~ e iq - y V {U), (101) 

then, the equations of motion become 

U^duiU^duAJ - R 3 q 2 A^ = , (102) 

R 3 q 2 A v - hj-^duiU^Au) - U^dUU^Au) = . (103) 
We obtain the solutions 

^=^M?rMw)' (i ° 4) 

^q-nq 1 ^ (R 3 yy* (2qR 3 l 2 \ 

where K 5 / 4 and Ki/ A are modified Bessel functions, and q = y/q 1 . The gauge fields satisfy 
the required boundary conditions vanishing rapidly for U < q 2 R 3 , due to the exponential 
behaviour of the Bessel functions. This implies that the perturbation will be more suppressed 
in the interior of the space defined by the 0123?7-coordinates as q increases. 
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DIS from scalar mesons 

In this section we show how to obtain the scalar mesons in the present model. Then, we 
obtain the interaction Lagrangian, and finally the structure functions for DIS from scalar 
mesons, using the method developed in the previous section. We begin with the metric 
induced on the D8-brane [5] 

ds 2 = (-J ^dy^dy" + (-J dZ ■ dZ , (106) 

where U —\Z\ and Z l (i — 1, • • • , 5) parameterize the 56789-space. Scalar mesons correspond 
to fluctuations of the D8-brane along the orthogonal directions. The dynamics of the D8- 
brane is described by the following action 

Sds = -Ha J d 9 ^-det(P[G] ab + 2na>F ab ) + J P[C&] A F A F . (107) 

The D8-brane tension is //§ = [(27r) 8 <7 s a' 9 / 2 ] -1 and P denotes the pullback of the background 
fields on the world-volume of the D8-brane. The EOM for scalar mesons arises from the 
transversal fluctuation 

x A = L + 2W<I> , (108) 
being $ the scalar fluctuation. From Eq.(107) one obtains the Lagrangian 



C = -/i 8V /|det(^ a6 )| v / det(5 afe + R^I 2 U^/ 2 g ac (2Tia') 2 d c ^d b ^) , (109) 

where all the indices denote directions on the D8-brane. Expanding the Lagrangian up to 
quadratic order we obtain 

'R^l 2 



£ = -^8Vl det (fl'' 



lab 



l + 2(^-J g a \na') 2 d a <$>d b <$> 



;no) 



Now, let us take spherical coordinates along the directions (Z l , ■ ■ •, Z 5 ), then use the metric 
(96), plug it into the quadratic Lagrangian, and obtain the following EOM 



d a (JdetgU 7 ^g ab d b ^)=0, (111) 



where is the ^-metric, which together with the U piece include directions (Z 1 , ■ ■ •, Z 5 ). 
Expanding the EOM around the interaction region one obtains 

R 3 d»d^ + U- 1/4 du(U 13/4 du^) + C/V^V^ = , (112) 

where Vj is the covariant derivative on S* 4 . Indices i, j, k, • ■ ■, denote coordinates on S* 4 , //, 
v, • • •, denote directions parallel to the D4-branes and m, n, ■ ■ •, correspond to 0, 1, 2, 3, U 
coordinates. The proposed Ansatz for the solution is 

$ = <t>(U)e ip - y Y e (S*), (113) 
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where Y t (S i ) are the spherical harmonics on S* 4 , 

V i V l Y e (S 4 ) = -£(£ + 3)Y\S*) . (114) 

Now, the EOM for becomes 

-P 2 R 3 <j) + U-^duiU 13 /^) - m\U<f> = , (115) 

with m 2 = £(£ + 3). For the initial hadronic state we can consider the leading behaviour of 
the solution for U ~ Ui nt where \R 3 Pf N / OUT \ <C U int . However, the intermediate state X 
must be represented by the full solution with P = P X - Thus, we obtain 

^in/out = Cl (Auy 9 ^ 2 e* p -yY(n) , (ne) 

$x = c x A- 1 ^ s V4^-9/8 J ^2 £ ^^^ l p x . ?/r(fi) (n7) 

where Jj, is the Bessel function of first kind, v = (1/4)^81 + 64£(£ + 3), being and q 
dimensionless constants. If we analyze the parity, we can classify the mesons as scalar and 
pseudoscalar mesons for even and odd values of £, respectively since [Y (fi)]p = (— 1) Y (fl). 

As in the D3-D7 case, the fields (113) are charged in the space 0123£7. All the fields with 
£ > have charge Q, since they transform with a phase under the U(l) symmetry group 
induced by transformations on S* 4 in the direction of the Killing vector. Thus, the quadratic 
Lagrangian from which we can derive the EOM (112) is 

C = -^a') 2 ^J\^\R z ' 2 U- 3 / 2 g ab {d a ^d b ^* + d h <$>d a <$>*) . (118) 

The U(l) current inserted at the boundary of the space spanned by the 0123?7-coordinates, 
produces a perturbation on the boundary conditions and excites a non-normalizable mode 
A m , which corresponds to the bulk extension of A^^- This field couples to the U(l) 
Noether's current in the bulk. The alluded perturbation propagates in the bulk inducing 
a metric fluctuation which is the product of A m (y,U) and a Killing vector Vi correspond- 
ing to the U(l) C 5*0(5) isometry of S* 4 . Therefore, we consider the fluctuation (26), 
Sg mi = A m (y,U)v l (VL). Using the eigenvalue equation v j djY(Q) = iQY(£l), we obtain the 
interaction Lagrangian 

C mt = iQ^{Tia') 2 ^\^\R 3 l 2 U- 3 / 2 A m ^d m ^ x - <$> x d m <$>) . (119) 

Now, let us compute explicitly C int from the coupling of the gauge field A rn to the Noether's 
current corresponding to the internal symmetry U(l) of the Lagrangian (118). Under an 
infinitesimal transformation with parameter -&/2 — constant we have 

j m = ifx 8 {na') 2 R 3/2 U- 3/2 (^d m ^ x - <^<9 m $) . (120) 
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Thus, if we define C int = Q\J\detg\A m j m , we obtain the same C int found in Eq.(119) from the 
metric fluctuation. Thus C int is indeed given by the coupling of a gauge field to a conserved 
current. 

In order to obtain the coupling A m j m , we use the expressions from (104) and (105), 

A. = ^^%T«»tfff). (12!) 



r(5/4) " 



A v = ~tT%QuA v , (122) 
qz 

and the current conservation 

d v f + U~ 13/4 du(U 13/4 j u ) = , (123) 

obtaining 

A m j m = Ap(j» - i^f 1 ) - i^U-^duiU^A.f) . (124) 
Now, we obtain the action 

S int = f°° d 9 xQ^tg~\A m3 m 

= |J"rf 9 x2/i 8 (W) 2 Qv^j^ 15/4 ^ 1/4 (^ + £)A^$x 

/ *V dQQ^/\d^\(A,U^Y)\T^ • (125) 

Using z = -R 3 / 2 /U 1 / 2 and 1/A — >■ oo, since A^ goes to zero very fast for U < U int , the result 
is 

(P + q,X\J»(0)\P,Q) =«s/ 3 (^ + £), (126) 

where q has been redefined, « 3 = 8/i 8 Q(7ra') 2 c { c\ [TX5/4)]-V /4 A , '~ 7/2 ar 1/4 (l -a;) 1 / 4 and 
the integral in [/ gives ^3, which using the s ~ g 2 (l/x — 1) approximation becomes 

J 3 = irf^ + u)q- v ~ WV/4 + ,/ 2(1 _ ^2 (127) 



4 V4 

In order to compute Im T M/3 from Eq.(23), we need to multiply Eq.(126) by its complex 
conjugate and proceed in an analogous way as before with the D3-D7 system. The distance 
between zeros of the Bessel functions in (117) is M n = nnA. The density of states is obtained 
in the same way as in the previous section. The final result of the tensor is 

ImT^ = 87r 5 Q 2 [r(9/4 + z/)] 2 [r(5/4)]-V 2 (a') 4 |Q| 2 |c x | 2 

A 2v -V 4 - 2 "(1 - x) v x^' 2 (V + £) (V + . (128) 



x 
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We obtain 

ft = 0, F 2 = E' Q 2 (^) (^)^ +7/2 (l " xY , (129) 
where E' = 47r 5 |cj| 2 |cx| 2 [r(9/4 + z/)] 2 [r(5/4)]~ 2 is a normalization dimensionless constant. 

DIS from vector mesons 

In this section we obtain vector mesons in the Sakai-Sugimoto model, compute the C int and 
derive the structure functions of DIS from polarized spin-one mesons in the context of the 
type IIA string dual description. In fact, we have obtained a full tower of vector mesons 
from the Dirac-Born-Infeld action of a D8-brane (Nf = 1) considering fluctuations along the 
brane, which has an additional quantum number £ with respect to the mesons calculated 
in [5]. We are specially interested in mesons with £ > as we need to study charged fields 
in order to be able to derive the interaction Lagrangian. So, for £ > 0, our mesons will 
transform with a phase under the U(l) symmetry induced by transformations on the S* 4 in 
the direction of a Killing vector. This can be satisfied by requiring the mesons not only to 
be dependent on the coordinates 0, 1, 2, 3, U, but they should also depend on the angular 
coordinates of the S* 4 . As a result, our mesons are charged and massive from the point of 
view of the five-dimensional space obtained after dimensional reduction on the S 4 . 

Besides, due to the fact that the relevant interaction region for the DIS corresponds to large 
values of the four-dimensional energy, we approximate the metric for that specific region. 
Consequently, we have calculated full analytic expressions for these mesons depending on 
the nine coordinates of the D8-brane. 

From the action (107) we derive the EOM, then we propose a quadratic Lagrangian from 
which we obtain exactly the same EOM, thus 

Cds = -^-det(g ab + 2na'F ab ) + ^£^ 8 P[C^] A F A F , (130) 

= C-DBI + C-CS ■ 

The Ansatz for the solution of vector mesons is 

B^C^(U)e ip - y Y e (S 4 ), P-C = 0, B u = 0, B t = , (131) 

where we have expanded it in spherical harmonics of S* 4 (there is a field for each £, as for 
the D3-D7 brane system). We have to determine 4>(U). The polarization vector is and 
P ■ ( = comes from d^B^ = 0. Taking into account the Ansatz above the EOM reads 

d a (J\d^\F ab ) = 0, (132) 
where a, b, c, d, ■ ■ ■ run over all the coordinates of the D8-brane. 
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Plugging the Ansatz in Eq.(132), we find that solutions with b = k and b = U identically 
satisfy the EOM from d^B^ = 0. If b = fj, we find 

R^U-S^^^PuBp-P^B^ + U'^^du^^duB^-U^^mlB, = 0, (133) 
which leads to an equation equal to Eq.(115) for (j>(U) while mj = £{£ + 3) comes from 

V 2 si Y e (S 4 ) = -£(£ + 3)Y e (S 4 ) . (134) 

The solution is 

m)=cU-^J v {^$^) ) (135) 

where J u is the Bessel function of first kind, c being a dimensionless constant, while v = 
(1/4)^81 + 64£{£ + 3). Again, as in the scalar case, for initial and final hadronic states we 
can use the leading behaviour for \P 2 \ <C q 2 , then 

B, in/ out = C»cA-\MJ)-^-V % e ip - y Y\S A ) , (136) 

being q a dimensionless constant. The intermediate state is given by the full solution 

B X » = C^cxA- 13 / 8 ^ 1 /^- 9 /^^^^)^^^ Y\S±) , (137) 

where cx is a dimensionless constant. The radius of the Bessel function 2s 1 / 2 i? 3 / 2 is much 
larger than Uq, being this the reason to consider the full solution for intermediate states. To 
study the parity transformations, using the relations [Y (fl)]p = (— l) e Y e (Q) and [CJp — 
(C°, — C) we can classify the mesons as vector mesons for even values of £ and axial vector 
mesons for odd values of £. 

From the expansion of B^ in spherical harmonics of S 4 , it can be seen that the gauge fields 
on the D8-brane correspond to charged fields in the space defined by 0123C7 coordinates in 
an analogous way as for the D3-D7-brane model. The quadratic Lagrangian which gives the 
EOM (133) for the solution (131) of vector mesons is 



C = -2^(na'Y^\detg\F ab F: b . (138) 

Considering a fluctuation of the metric of the form 5g m j = A m (y, U)vj(Q) with vWjY(Q) = 
iQY(Q) the interaction Lagrangian is: 

C = 2 l Q^{na'f^A^\A m [B x ^ m - B^)*] , (139) 

which can alternatively be obtained from the coupling of the gauge fiel d to th e Noether's 
current corresponding to the U{\) symmetry as before, being C int = Q\J\ det g\A m j m while 
the current is given by 

j m = 2i^a'f{B x ,F^ - B^FDl ■ (140) 
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Then, the coupling of A m to a conserved current j m leads to C int . 

In order to obtain the gauge field coupled to the conserved current, we use a similar 
procedure to that in the previous section, thus 

Sint = r d 9 xQj\d^\A mJ m 

JUo 



r d 9 x2fi 8 (7ia') 2 Q^\detg\R 3 U- 3 A^BB* x 

J Uo 

J d 4 ydnJ\d^\(A,U^Y)\%R3 , (141) 



'Uo 

gMjR 3/4 



q 2 



where = (^B and 

= 2(C • Cx) (P" + £) + (Cx • q)C ~ (C • q)Cx • (142) 

Considering z = R 3 ^ 2 /U 1 ^ 2 and 1/A — > oo since A^ goes to zero very fast for U < Vint-, 
the second term in Eq.(141) vanishes when evaluated at the integration limits. Thus by 
integration of the remaining term and comparing with Eq.(48), we obtain 

(P + q,X\J»(0)\P,Q) — K4I4N' 1 , (143) 

where k a = 8fisQ(7ra') 2 c i c x [T(5/A)]- 1 A u - 9 / 2 q 3 / i x 1 / 4 (l-x)- 1 / i and the result of the integral 
in U is I4, that in the approximation s~g 2 (l/x — 1) corresponds to 

h = 2- 2 r(v + 9/4)q-»- 13 / i x»/ 2+g /\l-xf 2 , (144) 

In order to compute Im T M/3 given in (23) we have to multiply Eq.(143) by its complex 
conjugate, add the radial excitations, and finally sum over the polarizations of the final 
hadronic states, obtaining 

A „l/2 



As 1 / 2 ^ 



7T fv^fv^l £ 



* t2 



E 



As 1 / 2 . 



i^vv)(/'" • I;) • kr -(C-ry)^ 



2(C* ' Cx) + ^) + (Cx • 9)^ - (C* • } • (145) 

Also, for the solution (131) it holds P-£ = 0, which implies that there are three polarizations 
A = 1, 2, 3. We use the normalization for the polarizations given by C M (Px, A) • (*(Px, A') = 
—M x 5x ; \', and we obtain 

J2(x»(Px,\)(h(Px,\) = M 2 xV ^ + P x ,P X p, 

A 

= M^ + (P M + ^)(P /9 + ^). 
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The rest of the computation is similar to the one done in the previous section for the D3-D7- 
brane system. We introduce the final form for the structure functions in the Sakai-Sugimoto 
model [5] which we have obtained: 





= E(x) 


——;{! -x-2xt- 4xH + 
Ylx 6 


4x 3 t + 8xH 2 ) , 


(146) 




= E(x) 


-^—(l-x + 12xt- UxH- 
ox 6 


- 12x 2 t 2 ) , 


(147) 




= E(x) 


— (1-x-xt), 




(148) 




= E{x) 


_L(l_ x _ a: 2 f ) > 




(149) 


h 


= E(x) 


24*3 (1 4 * + 8 ^)> 




(150) 




= E(x) 


Ux s( ! + 4 ^ 2 ^), 




(151) 


9i 


= E(x) 


-^-r(-3 + 3x + 4xt + 5x 2 t 


- QxH - 8xH 2 ) , 


(152) 


92 


= E(x) 


1 . (3 3a; Axt + 2x 2 t) 




(153) 



We can observe that they have the same form as the ones obtained with the D3-D7-model, 
however the factor E(x), which turns out to be common to all the structure functions, 
depends on the model and differs with respect to the previous model. E(x) is given by 

E & = ^oQ 2 (^) )> 11/2 (1 - -r 1 , (154) 
where E = 7r 5 4|cj| 2 |cx| 2 [r(z/ + 9/4)] 2 [r(5/4)]~ 2 is a dimensionless normalization constant. 

6 Results and Discussion 

In this section we discuss our results extensively. We have investigated deep inelastic scat- 
tering from scalar mesons and polarized vector mesons obtained from two holographic dual 
models with flavours in the fundamental representation of the gauge group. These mod- 
els are the D3-D7-brane model [1] and the D4-D8-D8-brane model [5], which are dual to 
the large N limit of SU (N) M = 2 supersymmetric Yang Mills theory and of SU (N) non- 
supersymmetric Yang-Mills theory within the same universality class as QCD 6 , respectively. 
Although we have only considered one flavour (Nf = 1), the method described above can be 
straightforwardly extended to the case of Nf > 1. 

6 There are certain differences with the large TV limit of QCD itself. For instance, the model in [5] has a 
global SO (5) symmetry which QCD does not have. 



30 



As discussed above, DIS amplitudes are obtained from the Lorentz contraction of a lep- 
tonic tensor with a hadronic tensor W^ v which characterizes the structure of the hadron. 
Since strongly coupled effects -which cannot be accounted for by using quantum field theory 
perturbative methods- contribute to this tensor, it is necessary to approach the problem in 
a completely different way. In fact, this is achieved through the implementation of a holo- 
graphic dual description, which makes it possible to investigate the strongly coupled regime 
of the field theory in terms of a string theory dual model. With this in mind, we have derived 
the hadronic tensor using two different holographic dual models, in the large N limit and for 
large values of the 't Hooft coupling. In this regime, as suggested by Polchinski and Strassler 
[10], the currents are scattered by the entire hadron rather than by any constituent parton. 
We shall comment more on this point. 

Let us briefly summarize a few properties satisfied by the holographic dual models we have 
considered. The interaction Lagrangian used in both models is invariant under parity and 
time reversal symmetries. Also we ought to mention that the existence of a conserved current 
in the holographic dual model implies that q^W^ = q^W^ = in the four- dimensional 
theory, where g M is the virtual photon four-momentum. In four dimensions this property 
comes from the conservation of the current coupled to the gauge field in the interaction 
Lagrangian. In order to obtain the interaction Lagrangian in the holographic dual models 
we have implemented a systematic method based upon the computation of the Noether's 
current of the DBI Lagrangian expanded up to quadratic order in the derivatives of the 
fields. By applying the method based on the Noether's current, we have found a C int with 
similar properties for both models, i.e. the D3-D7-brane model and the D4-D8-D8-brane 
model. This guarantees the existence of a conserved current and, thus q^W^" = q v W^ = 
is automatically satisfied. 

If we analyze the parity transformations in both models studied, we can classify the mesons 
in towers of scalar, pseudoscalar, vector and axial vector mesons. This classification depends 
on the values of £, which characterizes the spherical harmonics, Y e (Q), of the mesons wave 
function. It is interesting to stress that within the D4-D8-D8 model we have obtained 
analytic expressions for scalar as well as for vector mesons which are massive in the five- 
dimensional spacetime defined by the coordinates (y°, y 1 , y 2 , y 3 , U). Those expressions hold 
for large values of the energy of the dual four-dimensional quantum field theory, given by 
~ [J 3 / 4 //? 7 / 4 , which corresponds to U ^> Uo, where U$ = A 2 R 3 is the minimum of the 
coordinate U. The tower of vector fluctuations of the D8-brane that we have calculated has 
an additional quantum number, £, compared to the mesons considered in [5], that comes 
from the expansion in spherical harmonics on S A . 

We have obtained the DIS structure functions for the whole scalar and vector meson 
spectra using both holographic dual models, getting exact analytical expressions for these 
functions in all the cases. Particularly, for the vector mesons we have got the full tensor 
structure decomposition of the hadronic tensor for polarized targets of spin-1. In fact, we 
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have computed all the structure functions: F 1 , F 2 , b\, b 2 , 63, b 4 , g\ and g 2 , as functions of 
x and t. This problem had not been considered in the previous literature of vector mesons 
using the gauge/string duality, so this is a totally new approach to the problem which 
complements those based on the operator product expansion of electromagnetic currents 
and the parton model approach investigated by Hoodbhoy, Jaffe and Manohar [6]. The 
present calculation gives nontrivial results and predictions for the structure functions of 
spin-1 hadrons, particularly focussed on vector mesons. 

We find that Fi = for scalar mesons as expected. Recall that F± is proportional to 
the Casimir of the scattered hadron under the Lorentz group and therefore it vanishes for 
a scalar hadron as it would do for a scalar parton, in agreement with our result above. On 
the other hand, F 2 provides spin-independent information and does not vanish. Besides, for 
the scalar mesons, the dependence on momentum squared in the equation for F 2 , which is 
F 2 ~ g~ 2 ( T "p~ 1 ) ) corresponds to the twist r p = z/ + 2 in our notation (note that v is different in 
each model). For the vector mesons, we obtain finite results for all the structure functions, 
that we discuss below. 

Below we display our results about the structure functions in Figures (2-4), corresponding 
to the D3-D7-brane model (left) and D4-D8-D8-brane model (right). In all figures we have 
considered only the lowest values of the mass mi, which depends on the model. In these fig- 
ures only the dependence upon the Bjorken parameter is shown. In particular, the following 
functions are displayed 7 : 



F 2 x IO7 
F 2 x 10 6 / 



A 8 



2f m a 
V A 10 



A2 
q 2 
A2 

2 



(Fi, bj, 9i ) x 10 3 / 
(Fi, bj, 9i ) x 10 6 / 



A Q 2 
E Q 2 



A 8 
A 10 



A^ 
q 2 



for D3-D7 scalar mesons , 
for D4-D8-D8 scalar mesons , 
for D3-D7 vector mesons , 
for D4-D8-D8 vector mesons , 



(155) 



where i = l,2,j = 1, 2, 3, 4 and v = i+l for the D3-D7 model and v = (1/4)^81 + 64^ + 3) 
for the D4-D8-D8 model. In these figures we have neglected terms proportional to t with 
respect to the ones proportional to x n with n > 0, since the computation has been done 
within the regime where x ~ 1 (although horizonal axes include the full range < x < 1 we 
assume our calculations to hold for x ~ 1) and \t\ = (—P 2 /q 2 ) 1, since — P 2 q 2 (P is 
the ingoing meson four- momentum). 



7 In Figures (2-4) wc show the structure functions obtained in the previous two sections re-scaled as 
Eq.(155). Notice that we use the notation F\, F 2 , 61, b 2 , &3, &4, gi, g 2 for the corresponding re-scaled struc- 
ture functions. Also, the functions are multiplied by different factors (powers of 10 explicitly written in 
Eq.(155)) in each model in order to make the comparison between the two models easier. 
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For the scalar mesons, from Figure 2 it can be seen that the Bjorken parameter dependence 
of F 2 is similar in both models. This similarity between the structure functions in both 
models is extended to the vector mesons, as shown in Figures 3 and 4. The reason is that 
each structure function in both models shares a common factor which is model-independent 
and is determined by their corresponding DBI actions. However, each of them has a model- 
dependent factor, hence the differences between curves in both models. 

For the vector mesons, we have found the modified Callan-Gross relation as previously 
obtained in [10], namely 2^(2;) = F 2 (x), which holds for every value of the Bjorken param- 
eter in the present case. In reference [10] this modified Callan-Gross relation was obtained 
for the case of spin-| hadrons, corresponding to supergravity modes of the dilatino. In the 
present case we obtain the same relation for the supergravity modes of gauge fields which 
are fluctuations along the flat directions of the D7 and D8 brane world- volumes, depending 
on the model we consider. As pointed out in [10] the difference with respect to the usual 
Callan-Gross relation is an x-factor missing on the left-hand side of the relation. In the 
case of a scattered parton which carries a fraction x of the total four-momentum of the 
hadron, the operator product expansion leads to 2xF\(x) = F 2 (x). On the other hand, in 
the pure supergravity description the entire hadron is expected to be scattered due to the 
strong coupling limit in which the system is studied. We emphasize that within the super- 
gravity approximation our result 2Fi(x) = F 2 (x) is exact 8 . In addition, we also find the 
exact relation 2bi(x) = b 2 (x), which is the pure supergravity version of 2xbi(x) = b 2 (x), also 
obtained from the operator product expansion analysis of the time ordered product of two 
electromagnetic currents. The last relation comes from the fact that 61 and b 2 are obtained 
from the same operator tower as F 1 and F 2 , thus they obey the same scaling equations. 
Also notice that in the parton model F\ represents the probability of finding a quark in the 
hadron with momentum fraction x, while g\ is the difference between the probabilities of 
finding a quark in the hadron with momentum fraction x with spin parallel and antiparallel 
to the hadron [27]. Thus, in principle, g\ can be positive or negative. Our results agree with 
these physical interpretations. 

Also, we obtain the following relations: 26 3 = —64, 2g 2 = —gi, b 2 = 3F 2 and b\ = 3Fi, 
which we regard as predictions of both holographic dual models in the t — > limit. In 
particular the relation bi = 3Fi is in agreement with the expectations for a spin-one hadron 
since b\ ~ O(Fi), as in the case of the p- meson, consisting of relativistic spin-half constituents 
[6]. We would like to point out the fact that, on the one hand, within the supergravity 
approximation the structure functions do not reveal the partonic structure of the hadron 
(since there is a missing x-factor in the Callan-Gross relations), and on the other hand from 
bi ~ O(Fi) one may infer that the meson is composed by relativistic spin-half constituents. 
A possible explanation for this behaviour is that, as suggested in [10], at large 't Hooft 

8 For this relation as well as for the ones that are considered below, terms proportional to t have been 
neglected. 
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Figure 2: Structure function F 2 {x) for scalar mesons, scaled as explained in the text, as a function 
of the Bjorken parameter. Left and right figures correspond to the D3-D7-brane and the D4-D8-D8- 
brane models, respectively. Different curves stand for different values of t > 0, we display the lowest 
values of it. The curve sizes increase with t. Notice that t has been neglected. 



coupling all partons are wee since parton evolution is so rapid that the probability of finding 
a parton with a substantial fraction of the incoming momentum of the initial hadronic state 
is vanishingly small. So, this explains the modified Callan-Gross relations. Now, at finite TV 
there is an argument to think of the parent hadron as surrounded by a cloud of other hadrons, 
and therefore the lepton should be scattered by one of the hadrons in the cloud. So, perhaps 
this is what b\ = 3F± indicates in the strongly coupled and large Bjorken parameter regimes. 
We do not have a completely satisfactory answer for this apparent contradiction between 
no partons seen from one type of relations and the existence of relativistic components of 
spin-1/2, which is inferred from the other relation. Also, we should notice that from OPE 
analysis [6] it is expected that F 1 , F 2 , &i, b 2 and g\ are twist-two structure functions, while 
the rest should not receive contributions at leading twist. Thus, this seems to be a sort of 
discrepancy with respect to our results where all the structure functions are of comparable 
order of magnitude. Moreover, in each model all the structure functions have a similar 
dependence of momentum squared, i.e. g -2 ^ -1 ), where the twist r p is u+ 1 in our notation 
(note that v is different in each model). It would be very interesting to understand these 
differences, in order to establish whether the considered holographic dual models are good 
enough as to describe the internal structure of realistic vector mesons, including both the 
polarized and the unpolarized contributions. 

Extensions of this work can be done to more general holographic dual models. In partic- 
ular, it should be straightforward to apply these ideas to the D4-D6-D6 model [4], as well as 
to extend the models considered here to the non-Abelian cases, i.e. Nf > 1. Also it would 
be interesting to apply this approach to the case of flavoured versions of the Klebanov- 
Strassler solution [29], as in the cases of [30] and [31]. Likely, flavoured versions of the 
non-supersymmetric models presented in [32] and [33] might be developed and DIS studied 
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Figure 3: Structure functions F\[x), i*2(x), b\{x) and b2(x) for vector mesons as a function of 
the Bjorken parameter scaled as explained in the text. Left and right figures correspond to the 
D3-D7-brane and the D4-D8-D8-brane models, respectively. Different curves stand for different 
values of £ > 0, we display the lowest values of it. The curve sizes increase with I. Notice that t 
has been neglected. 
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Figure 4: Structure function 63 (x), 64 (x), gi(x) and gi (x) for vector mesons as a function of the 
Bjorken parameter scaled as explained in the text. Left and right figures correspond to the D3-D7- 
brane and the D4-D8-D8-brane models, respectively. Different curves stand for different values of 
£ > 0, we display the lowest values of it. The curve sizes increase with £. Notice that t has been 
neglected. 
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in these cases. Similarly, these ideas could likely be also extended to the Maldacena- Nunez 
solution [34] with the addition of matter fields in the fundamental representation [35]. For 
all these models it would be very interesting to unveil how the DIS structure of spin-one 
mesons is manifested. 
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